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> \ 1 INTRODUCTION 

• i-H , 

It is well known that a stratification in the chemical com- 
■ position of stellar models directly influences the properties 
of gravity modes. The signatures of chemical stratifications 
have been extensively investigated theo retically, and ob- 
served, in pulsating white dwarfs (see e.g. iKawalcr 19951, for 
a review). The influence of chemical composition gradients 
on g-modes in main sequence stars has been partly addressed 
and sugges ted in the works bylBerthomieu fc Provosti (1988) 
and lDziembowski et alj (1 19951). Inspired by these works and 
following the approac h of Berthomieu fc Provosti (| 19881 ) and 
iBrassard et al.l (I1992T I we study the properties of high-order, 
low-degree gravity modes in main sequence stellar models. 

High-order gravity modes are observed in two classes 
of main-sequence stars: 7 Doradus and Slowly Pulsating 
B stars. The former ar e main-sequence s tars with masses 
around 1.5 M (see e.g. IGuzik et al.ll2000h that show both 
photometric and line profile variations. Their spectral class 
is A7-F5 and their effective temperature is betwe en 7200- 
7700 K on the ZAMS and 6900-7500 K above it (|Handlerl 
1999). 7 Dor stars are multiperiodic oscillators with periods 
between 8 hours and 3 days (high-order g-modes). The mod- 
ulation of the radiative flux by convection at the base of the 
convective envelope was proposed as the excitation mecha- 



ABSTRACT 

In main sequence stars the periods of high-order gravity modes are sensitive probes 
of stellar cores and, in particular, of the chemical composition gradient that develops 
near the outer edge of the convective core. We present an analytical approximation 
of high-order g modes that takes into account the effect of the ^ gradient near the 
core. We show that in main-sequence models, similarly to the case of white dwarfs, 
the periods of high-order gravity modes are accurately described by a uniform period 
spacing superposed to an oscillatory component. The periodicity and amplitude of such 
component are related, respectively, to the location and sharpness of the /1 gradient. 

We investigate the properties of high-order gravity modes for stellar models in a 
mass domain range between 1 and 10 M , and the effects of the stellar mass, evo- 
lutionary state, and extra-mixing processes on period spacing features. In particular, 
we show that for models of a typical SPB star, a chemical mixing that could likely 
be induced by the slow rotation observed in these stars, is able to significantly change 
the g-mode spectra of the equilibrium model. Prospects and challenges for the aster- 
oseismology of 7 Doradus and SPB stars are also discussed. 
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nism for such stars (see e.g. IGuzik et al. l l2000l : lDupret~ et al.l 
|2004) . 

Slowly Pulsating B stars (SPB) are multiperiodic main- 
sequence stars wit h masses from a bout 3 to 8 Mq and spec- 
tral type B3-B8 jWaelkensl Il99ll ). High order g-modes of 
periods typically between 1 and 3 days are found to be ex- 
cited by the K-mechanism acting in the region of the metal 
opacity bump lo cated at T ~ 2 x 10 5 K in the stellar 
interior (see e.g. Dzicm bowski et al] Il993]). Recent obser- 



vations (see Jerzvkiewicz et al. 20051 : Handler et al.1 120061 : 



IChapellier et all |2006h and theoretical instability analysis 



l Pamvatnvkhlll999l : iMiglio et al.l 120071 ) also suggest high- 
order g-modes being excited in a large fraction of the more 
massive j3 Cephei pulsators: the seismic modelling of these 
hybrid pulsators looks very promising as it would benefit 
from the information on the internal structure carried by 
both low order p and g modes (/3 Cephei type oscillation 
modes) and high order g modes (SPB-type pulsation). 

The seismic modelling of 7 Doradus and SPB stars 
is a formidable task to undertake. The frequencies of 
high-order g modes are in fact closely spaced and scan 
be severely perturbed by t he effects of rotati on (see e.g. 
iDintrans fc Rieutordl |2000| : ISuarez et all [2005). Nonethe- 
less, the high scientific interest of these classes of pul- 
sators has driven efforts in both the observational and the- 
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oretical domain. Besides systematic photometric and spec- 
troscopic ground-ba sed surveys carried ou t on 7 Dor (see 
Mathias et all [2004) and SPB stars (see iDe Cat fc Aertsl 



0.8 



2002i ) , the long and unint errupted photomet r ic observations 



plann ed with COROT l|Baglin et all 120061 ; iMathias etail 
120061 ) will allow to significantly increase the number and 

accuracy of the observed frequencies. 

On the theoretical side, as suggested by ISuarez et al.l 

(2005) in the case of 7 Dor stars, a seismic analysis becomes 
feasible for slowly rotating targets. In th ese favorable c ases 
the first-order asymptotic approximation (|Tassoul|ll98(]| ) can 
be used as a tool to d erive the buoyancy radius of the star 
(see iMova et alj|2005h from the observed frequencies. Nev- 
ertheless, the g-mode spectra of these stars contain much 
more information on the internal structure of the star. In 
this paper we describe in detail the information content car- 
ried by the periods of high-order g modes, and show that 
the effect of chemical composition gradients can be easily 
in cluded a s a re finement of the asymptotic approximation 
of lTassoull (|l980h . 

After an introduction to the properties of gravity modes 
in main-sequence stars (Sec[2]), we present in Sec.|3]the ana- 
lytical approximation of high-order g-mode frequencies that 
will be used in the subsequent sections. In Sec.|3]we describe 
the properties of numerically computed g-mode frequencies 
in main-sequence stars in the mass domain 1-10 Mq. The 
effect of adding extra-mixing at the outer edge of the con- 
vective core (rotationally induced turbulence, overshooting, 
diffusion) is investigated in Sec. 14.21 In Sec. [5] we estimate 
how of the effects of rotation and of current observational 
limitations affect asteroseismology of main sequence high- 
order g modes pulsators. A summary is finally given in Sec. 

El 



2 THE PROPERTIES OF TRAPPED g-MODES 

As it is well known, the period spectrum of gravity modes is 
determined by the spatial distribution of the Brunt- Vaisala 
frequency (TV) which is defined as: 



TV 



1 dp 1 dp 
Fip dr p dr 



(1) 



TV can be approximated, assuming the ideal gas law for a 
fully-ionized g £LS j cLS '. 



TV 



9 2 P 



(Vad - V + V„ 



where 



V = 



dlnT 
d lnp ' 



V ad = 



dlnT 
d lnp 



and Vu = 



d ln/x 
d lnp 



(2) 



(3) 



The term gives the explicit contribution of a change of 
chemical composition to T V. The first or der asymptotic ap- 
proximation developed bv lTassoull (|l980h shows that, in the 
case of a model that consists of an inner convective core and 
an ou ter radiative envelope (we refer to the work bv lTassoull 
1 19801 . for a complete analysis of other possible cases), the 
periods of low-degree, high-order g modes are given by: 



Lf 1 

Jx 



'-dx 



(2k + n e ) 



(4) 




m/M 



Figure 1. Upper panel: Hydrogen abundance in the core of 6 
Mq models on the main sequence. X c ~ 0.5 (dotted line) and 
at X c ~ 0.3 (dashed line). The convective core recedes during 
the evolution and leaves behind a chemical composition gradient. 
Lower panel: Buoyancy frequency N as a function of the normal- 
ized mass. 



where L = [£(£ + l)] 1 ^ 2 (with I the mode degree), n e the 
effective poly tropic index of the surface layer, x the nor- 
malized radius and xq corresponds to the boundary of the 
convective core. In order to avoid confusion with n e , the 
radial order of g modes is represented by k. 

Following Eq. f?] the periods are asymptotically equally 
spaced in k and the spacing decreases with increasing L. It 
is therefore natural to introduce, in analogy to the large fre- 
quency separation of p modes, the period spacing of gravity 
modes, defined as: 



AP 



fe+i 



Pk 



(5) 



In the following sections we will show that deviations from a 
constant AP contain information on the chemical composi- 
tion gradient left by a convective core evolving on the main 
sequence. 

We consider as a first example two models of a 6 
Mq star evolving on the main-sequence. The behaviour of 
TV and of the chemical composition profile is represented in 
Fig. [Tj The convective core is fully mixed and, therefore, the 
composition is uniform (V M = 0). However, in stars in this 
mass range, the convective core shrinks during the evolution, 
leaving behind a steep gradient in the hydrogen abundance 
X. This causes a sharp peak in and in TV: does this fea- 
ture leave a clear signature in the pr operties of g-modes? 

This question was addressed by iBrassard et al l (Il99ll ) 
while studying the seismic effects of compositional layering 
in white-dwarfs. The authors found that a sharp feature in 
the buoyancy frequency could lead to a resonance condition 
that may trap modes in different regions of the model. 

A first indicator of such a trapping is the behaviour of 
(x) defined by: 



(x) 



J x\5r\ 2 dm 
J Q \Sr\ 2 dm 



(6) 
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Figure 2. (a;) for 6 Mq models on the main sequence at two 
different points in the evolution: X c ~ 0.5 (dotted lines) and at 
Xc — 0.3 (dashed line). Modes of different radial order are peri- 
odically trapped in the region of chemical composition gradient. 
Trapped modes are marked with open symbols. 



where x = r/R and Sr is the total displacement vector. 
As shown in Fig. [5] modes of different radial order k are 
periodically confined closer to the center of the star. 

In Fig. [3] we show the behaviour of the eigenfunctions 
for modes of radial orders around a trapped mode: the partly 
trapped mode has, compared to "neighbour" modes, a larger 
amplitude in the region of mean molecular weight gradient. 

In white dwarfs it has been theoreticall y predicted and 
obser ved (see, for instance, the recent work of Metcalfe ct al. 
2003) that the period spacing AP(fc) = Pfc+i — Pk is not 
constant, contrary to what is predicted by the first order 
asymptotic approximation of gravity modes. This has been 
interpreted as the signature of chemical composition gradi- 
ents in the envelope and in the core of the star. In analogy 
with the case of white dwarfs, in models with a convective 
core, we expect the formation of a nonuniform period dis- 
tribution; this is in fact the case as is presented in Fig. [4] In 
that figure we plot the period sp acing derived by using the 
adiabatic oscillation code LOSC l|Scuflaire et al. l2007al ) for 
models of 6 Mq at two different stages in the main sequence 
evolution. The period spacing presents clear deviations from 
the uniformity that would be expected in a model without 
sharp variations in N. How these deviations are related to 
the characteristics of the chemical composition gradients will 
be studied in the following sections. 



3 APPROXIMATE ANALYTICAL 

EXPRESSION OF g-MODES PERIOD 
SPACING 

In this section we derive two approximate expressions that 
relate deviations from a uniform period distribution to the 
characteristics of the /^-gradient region. These simplified ex- 
pressions could also represent a useful tool to give a direct in- 
terpretation of an observed period spectrum. Though a first 
description of these approximated expressions was outlined 
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Figure 3. We consider a 6 Mq model with X c ~ 0.5. Upper 
panel: N against fractional radius. Central panel: Horizontal com- 
ponent of the displacement for three £ = 1 gravity modes of radial 
order 18 (dotted line), 20 (continuous line) and 22 (dashed line). 
The eigenfunction corresponding to k = 20 is partly trapped in 
the region of mean molecular weight where the sharp variation of 
JV is located. Lower panel: As in the central panel, this time the 
radial displacement is shown. The eigenfunctions corresponding 
to different modes are normalized to have the same total pulsation 
energy. 



in lMigliol (|2006h and lMiglio et al.1 (|2006t ). we here present a 
more detailed analysis. 

We recall that deviations from the asymptotic ex- 
pressions of the frequencies of high-order pressure modes 
have been widely studied in the context of helioseis- 
mology. The oscillatory features in the oscillation spec- 
trum of solar oscillation modes allowed modes to infer 
the properties of localized variations of the solar struc- 
ture, e.g. at the base of the convective envelope and 



in the 


second helium ionization region (see e.g. Cough 


1990 


: 


Christensen-Dalsgaard et al. 


199ll: Monteiro et al. 


1994 


Basu & Antial Il995: Montciro & Thompson 2005: 


Houdek & Cough 20071). 



3.1 Variational principle 

A first and simple approach to the problem is to make use of 
the variational princip le for adiabatic stellar oscillations (see 
e.g. lUnno et al.lll989l ). The effect of a sharp feature in the 
model (a chemical composition gradient, for instance) can 
be estimated from the periodic signature in SP, defined as 
the difference between the periods of the star showing such 
a sharp variation and the periods of an otherwise fictitious 
smooth model. 

We consider a model with a radiative envelope and a 
convective core whose boundary is located at a normalized 
radius xq. N- and N+ are the values of the Brunt- Vaisala, 
frequency at the outer and inner border of the /^-gradient 

region. We define a — (^W^ with iV-|- iV_. Then a = 1 
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Figure 4. Period spacing for the same 6 Mq models as in Fig. [2] 
The periods of the components (in terms of k) are approximately 
7 and 3. Horizontal dotted lines represent constant period spacing 
predicted by the asymptotic approximation (Eq. [4J. 



describes the smooth model and a — > a sharp discontinuity 
in N. 

To obtain a first estimate of 6P, we adopt (following 
the approach by iMontgomerv et alllioolh the Cowling ap- 
proximation, that reduces the differential equations of stellar 
adiabatic oscillations to a system of the second order. Fur- 
thermore, since we deal with high-order gravity modes, the 
eigenfunctions a re well descr ibed by their JWKB approxi- 
mation (see e.g. Gough 1993). We can therefore express 5P 
as: 



8P k 



2IIo 



l ° fSN\ (LP k 7T 

cos 



, , cos ( + — ) dIl T 

V N J \tyU x 2/ 



(7) 



where L = [£(£+ 1)] 1/2 , the local buoyancy radius is defined 

\N\ 



nr 1 



-dx' , 



and the total buoyancy radius as 



Hn 1 



\N\ 



dx' 



(8) 



(9) 



(10) 



We model the sharp feature in located at x = x M 



The buoyancy radius of the discontinuity is then: 

n - = r ^ . 



as: 

5N _ I -a 2 
~W ~ a 2 



(11) 



where H(x) is the step function (see left panel of Fig. [5]). 

Retaining only periodic terms in 5P and integrating by 
parts we obtain: 



X]D IIo 1-a 2 / LP k n 

8P k oc — — cos — — + - 

La 2 \ ttH.i 2 



(12) 



N.(a.) 



N_(a,) 




N.(a 2 ) ; 



N_(o,) = 



Figure 5. Left: We model the sharp variation of N at x = x M by 
means of Eq. JTTJ for different values of a: a| = 0.5 (continuous 
line), a| = 0.3 (dashed line). Right: A smoother variation of N is 
modeled following Eq. 1151 . 



For small SP we can substitute the asymptot i c app roxima- 
tion for g-modes periods derived by iTassoull l|l980t ) in the 
expression above: 



P k =vr 



2 Ho 
L 



'2k + 4>') 



where <f>' is a phase constant that depen ds on the bou ndary 
conditions of the propagation cavity (see lTas soul 198fj), and 
find 



^oc^I^cos( 2 A + "V + f 
La 2 - \ IF ll u 2 



(13) 



From this simple approach we derive that the signature of a 
sharp feature in the Brunt- Vaisala, frequency is a sinusoidal 
component in the periods of oscillations, and therefore in 
the period spacing, with a periodicity in terms of the radial 
order k given by 



Afc 



Ik 

n 



(14) 



The amplitude of this sinusoidal component is proportional 
to the sharpness of the variation in iV and does not depend 
on the order of the mode k. 

Such a simple approach allows us to easily test the ef- 
fect of having a less sharp "glitch" in the Brunt- Vaisala 
frequency. We model 5N (Fig. [5] right panel) as a ramp 
function instead of a step function: 



5N 

~w 



a 2 (xy, - x) 

• 2 Xn — XQ 



(15) 



In this case integration by parts leads to a sinusoidal 
component in SP k whose amplitude is modulated by a factor 
1/Pfc and therefore decreases with increasing k, i.e. 



X]D i noi-Q 2 i / no, , n„ 

oP k OC — ^ , COS I 2lV -^—k + 7T- 



P k L 



nr, 



n. 



+ 



.(16) 



The information contained in the amplitude of the sinusoidal 
component, as will be presented in Section[4]2] is potentially 
very interesting. It reflects the different characteristics of the 
chemical composition gradient resulting, for example, from 
a different treatment of the mixing process in convective 
cores, from considering microscopic diffusion or rotationally 
induced mixing in the models. 

Eq. (|13p was derived by means of a first order perturba- 
tion of the periods neglecting changes in the eigenfunctions. 
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This approximation is valid in the case of small variations 
relative to a smooth model, therefore it becomes question- 
able as the change of N at the edge of the convective core 
becomes large. A more accurate approximation is presented 
in the following section. 

3.2 Considering the effects of the p gradient on 
the eigenfunctions 

We present in this Section a description of mode trapping 
considering the change in the eigenfunctions due to a sharp 
feature in N. As a second step we derive the effects on the 
peri ods of g-modes. 

iBrassard et~aH l|l992l ) studied the problem of mode 
trapping in ^-gradient regions inside white dwarfs. In this 
section we proceed as IBrassard et all (1 19921) , apply ing the 
asymptotic theory as developed in ITassoull l|l98ol ) to the 
typical structure of an intermediate mass star on the main 
seq uence. 

ITassoull (|l980h . assuming the Cowling approximation, 
provided asymptotic solutions for the propagation of high- 
order g-modes in convective and radiative regions, located 
in different parts of the star. In order to generalize the ex- 
pression for these solutions, she introduced two functions Si 
and S2 related to the radial displacement (£ r ) and pressure 
perturbation (p') as follows: 



2f 



and 



-1/2 -2i 
P x 



-1/4 



<?£(£+ l)p' /p = p 



Si 



-V2UI1/4 



s 2 



(17) 



(18) 



where x is the normalized radius, a the angular frequency 
of the oscillation, and 



'+ l)iV 2 



(19) 



The expressions for Si, & for different propagation re- 
gions are given in the equations [T79] to [T97]j. 

The only diff erence from the derivation of 
IBrassard et~afl ljl992h . who assumed an entirely radia- 
tive model, is that here we consider a model that consists 
of a convective core and a radiative envelope. The solutions 
should then be described by [T80] close to the center, [T96] 
in the convective core (xq), [T97] in the radiative region 
with [T82] close to the surface (where the structure of the 
surface layers of the model is described by an effective 
polytropic index n e ). 

Now, as explained in Fig. [6] we assume that at x = x^ 
in the radiative zone there is a sharp variation of N due to 
a p gradient and, as in the previous section, we model it as 
a discontinuity weighted by a, where a = (JV+ /iV_) (see 
Eq-EJl. 

We define A = cr -1 , 

rl \N\ 



v (x) 



l -dx' = L (iln 



n: 



and vi(x) = L I ^j-dx' = LIZ,, 1 

Jx x 



1 Henceforth on [Ti] indicates equation number i in the paper by 
ITassoull (|1980T) . 



Asymp. (S^SJ, 







Figure 6. A schematic view of the simplified model we consider. 
The radiative region outside the convective core (x > xq) is di- 
vided in two zones: one below and the other above the disconti- 
nuity in N located at x = x^. In each of these regions we consider 
the asymptotic expressions for the eigenfunctions that are then 
joined continuously at x = x^. 



For large values of Auo and Xvi, we can write the eigen- 
functions in the radiative region above the discontinuity at 

X ^ cLS '. 

Sia oc k sin (A) and fea oc — k cos (A) , (20) 
and in the Vp region below the discontinuity we have: 
Sib oc k\ cos (B) and 52b oc —ki sin (B) , (21) 
where 

A = \v {x)-n e ^- - A , B = \vi(x) - ^ , (22) 



and ko and ki are arbitrary constants. Note that A and B 
are functions of x. 

The eigenfrequencies are now obtained by matching 
continuously the individual solutions in their common do- 
main of validity. In particular, imposing the continuity of p' 
and £ r at the location of the discontinuity in N w e obtain 
the following conditions (as in lBrassard et af1ll992h : 



aS 1 , 
= 5 2 " , 



(23) 
(24) 



where ft, 2 is evaluated above (<S^" 2 ) an< l below (S 1 . 



This finally leads to a condition on the eigenfrequencies 



1/A: 
cos (A + B) = 



1 + a- 



■ cos (A - B) 



(25) 



With A and B fEq. I22|l evaluated at x = x^, this condition 
can also be explicitly written as: 



1-Q J 

1 + a 2 



cos XL 



aA_ 
no 


n e 7r 




2 




f) 


" 1 


2 " 




n e 7r 


ThT~ 






~2~ 



(26) 



3.2.1 Further approximations 

A first extreme case for Eq. (|26[) is that corresponding to 
a — 1, i.e. no discontinuities in the Brunt- Vaisala, frequency. 
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In this case Eq. ((25} immediately leads to the condition 
cos (A + B) =0 and therefore to the uniformly spaced pe- 
riod spectrum predicted by Tassoul's first order approxima- 
tion: 

Pk = n 2 ^(2k + n e ) (27) 

Another extreme situation is a — » 0, in this case N is so 
large in the ^-gradient region that all the modes are trapped 
there; the periods of "perfectly trapped" modes are then: 

P„ = (n+i)2^ 2 ^ (28) 

where n = (1, 2, 3, . . .). 

The interval (in terms of radial order fe) Afe between 
two consecutive trapped modes (An = 1) c an be obtained 
comb ining Eq. (|27p and Eq. (|28p (see also iBrassard et al.l 
1992), and is roughly given by: 

Afc~5& (29) 
Ho 

which corresponds to Eq. (|14p . 

We choose the 6 Mq models considered in Sec. [2] to 
compare the g-mode period spacings predicted by equations 
(|14[) and (|26|) . with the results obtained from the frequencies 
computed with an adiabatic oscillation code. 

Equation (I14|l relates the period of the oscillatory com- 
ponent in the period spacing AP to the location of the 
sharp variation in TV. In Fig.|4]the periods (in terms of fc) of 
the components are approximately 7 and 3 for models with 
X c = 0.5 and 0.3. Following Eq. (|14[) these periods should 
correspond to a location of the discontinuity (expressed as 
rio/TIu ~ fc" 1 ) of 0.14 and 0.3: as shown in Fig. these 
estimates describe very accurately the locations of the sharp 
variation of N in the models. 

Numerical solutions of Eq. (126b . found usi ng a 
bracketing-bisection method (see iPress et al.l [l992), are 
shown in Fig. [8] As is clearly visible when comparing Figs. 
[8] and [4l we find that the solutions of Eq. (|26[1 better match 



the oscillatory behaviour of the period spacing than the si- 
nusoids of Eq. (fT3)l . 

In the following section we extend to a wider range 
of main-sequence models the analysis presented for a 6 
M Q model. 



4 APPLICATION TO STELLAR MODELS 

The occurrence of a sharp chemical composition gradient in 
the central region of a star is determined by the appearance 
of convection in the core and by the displacement of the con- 
vective core boundary during the main sequence. For a given 
chemical composition, and if no non-standard transport pro- 
cess is included in the modelling, the transition from radia- 
tive to convective energy transport, as well as the shape of 
the ii gradients in the central stellar region, are determined 
by the mass of the model. On the other hand, additional 
mixing processes may alter the evolution of the convective 
core and the detailed properties of the chemical composition 
profile. 

As shown in the previous section, the features of pe- 
riodic signals in the period spacing of high order g-modes 
can provide very important information on the size of the 
convective core and on the mixing-processes able to change 
the n gradients generated during the evolution. 

In this section, we present a survey of the properties of 
adiabatic I = 1 high order g-modes in main-sequence stars 
with masses from 1 to 10 Mq, and for four different evo- 
lutionary stages: those corresponding to a central hydrogen 
mass fraction X c of 0.7, 0.5, 0.3 and 0.1. All these models 
were computed with the same initial chemical composition 
(Xq, Zo) = (0.70,0.02). The adiabatic oscillation fr equencies 
were computed with LOSC IjScuflaire et al.ll2007al ). 

We first study how the properties of high order gravity 
modes depend on the mass and the evolutionary stage of the 
model (Sec. 14. 1) . In a second step we evaluate the effects of 
the inclusion of extra-mixing such as overshooting, diffusion 
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Figure 9. Fractional mass of the convective core as a function 
of central hydrogen abundance for models computed with (right 
panel) and without (left panel) overshooting, for masses between 
1.0 and 2 Mq. (3 is the overshooting parameter as defined in 

Sect, rami 



and turbulent mixing (Sec. 14. 2p . The behaviour of modes 
with different I will be briefly addressed in Sec. 14.31 



4.1 Convective core evolution: stellar mass 
dependence 

In our analysis of the signatures of the fi gradients on g- 
modes, we consider three stellar mass domains: i) M < 
Mlcc, with Mlcc being the minimum mass required to keep 
a convective core during the main sequence; ii) Ml C c^ M ^ 
Mg C , for which the mass of the convective core increases dur- 
ing part of the main sequence; and Hi) M > M gc for which 
the convective core recedes as the star evolves. The situa- 
tions described here above are presented in Fig. |5J where the 
size in mass of the convective core is shown as a function of 
the central hydrogen abundance (and therefore of the age) 
of the star. The exact values of Ml cc and M > M gc depend 
on the chemical composition and, as we shall see below, on 
extra-mixing processes. 



4-1.1 Models with a radiative core 

As a first example we consider the evolution of the period 
spacing on the main sequence in models without a convec- 
tive core, e.g. in a 1 Mq star. The behavior is substantially 
different from higher mass models: AP, as shown in Fig. 
I lOt considerably decreases during the main sequence: this 
can easily be understood recalling the first order asymp- 
totic expression for the mean period spacing. The increase 
of N near the centre of the star, due to the mean molecular 
weight gradient developing in a radiative region (see upper 
panel of Fig. 1 1 Of) . has a larger and larger contribution to 
J §rdx' , leading to a significant reduction of the mean pe- 
riod spacing. The increase of N near the centre is however 
not sufficient to produce any periodic component of appre- 
ciable amplitude in the period spacing (see lower panel of 
Fig. [10). 
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Figure 10. Upper left panel: N as a function of the normalized ra- 
dius in 1 Mq models with decreasing central hydrogen abundance. 
Upper right panel: Hydrogen abundance profile versus r / R. Lower 
panel: g-modcs period spacing as a function of the radial order k. 
The asymptotic value of AP (predicted by Eq. @J| is represented, 
for each model, by dotted lines. 



4-1.2 Models with a growing convective core on the 
main-sequence. 

In models with masses between Mlcc and M gc , the con- 
tribution of nuclear burning through CNO cycle becomes 
more and more i mportant as the star evolves on the main 
sequence (see e.g. iGabriel fc Noelslll977l : ICrowe fc Mitalasl 
ll982l : |Popielski fc Dziembowskill2005t) . The ratio L/m in the 
nuclear burning region becomes large enough to alter the be- 
haviour of V ra d: the latter increases and so does the size of 
the convective core. 

A growing convective core generates a discontinuity in 
the chemical composition at its boundary (see Fig. Hip , and 
may lead to an inconsistency in the way the convective 
boundary is defined. The situation is illustrated in Fig. 1121 
the discontinuous hydrogen profile forces the radiative gra- 
dient to be discontinuous and to increase outside the region 
that is fully mixed by convection, and therefore, this region 
should be convective as well! If this is the case, then we have 
a contradictory situation: if we allow this region to have the 
same chemical composition as the core, then V ra a decreases 
and the region becomes radiative again. The question of the 
semi-convection onset in models w ith masses in th e range 
1.1- 1.6 was already address ed by IGabriel fc Noels! (|l977l ) 
and[c rowc fc Mitalasl (|l982h quite some time ago. Never- 
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Figure 11. The discontinuous chemical composition profile gen- 
erated by a growing convective core in 1.3 Mg (see Fig. 1126 when 
no extra-mixing outside the core is allowed. 
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Figure 12. Upper panel: Radiative gradient in the inner regions 
of 1.3 Mq models at different stages on the main sequence. The 
adiabatic gradient is represented with a short-dashed line. Middle 
panel: Ratio L(r)/m(r). During the evolution the large increase 
of L/m at the former border of the convective core dominates the 
behaviour of V ra d- The vertical dashed line denotes the border of 
the convective core in the ZAMS model. Lower panel: Behaviour 
of k for the same models as in the other panels. 



theless, what happens in this so-called "semi-convective" 
region is still a matter of debate. Some mixing is likely to 
take place, so that the composition gradients are adapted to 
obtain V ra d = V a( j in the semiconvective region. 

Even if no specific mixing is added in the semi- 
convective region, the \i gradient at the boundary of the 
convective core is very sensitive to the details of the nu- 
merical algorithm used in describing the core evolution. In 
fact, a strict discontinuity in chemical composition is only 
obtained if the border of the convective region is treated 



with a double mesh point (Fig. [TT] for instance). This "un- 
physical" framework leads to a problem when computing 
the Brunt- Vaisala frequency. The numerical difficulty can 
however be avoided keeping a quasi-discontinuous chemical 
composition, with a sharp change of X in an extremely nar- 
row region (So = — xq) outside the convective core (xo)- 
From Eqs. (|10[l and (|13[) it is evident that the signal in the 
period spacing will then have an almost infinite period. 

Of course, any treatment of the semi-convective re- 
gion should destroy the discontinuity leading to a wider 
pi-gradient region. The chemical composition discontinuity 
may also be removed by a sort of "numerical diffusion" that 
appears when the grid of mesh points (necessarily finite) in 
the modelling does not follow the con vection limits. That is 
the case of the evolution code (CLES. IScuflaire et al . 2007b) 
used to compute most of the stellar models presented in this 
paper. In these models, the region where the discontinuity 
would be located is assumed to have an intermediate chem- 
ical composition between the one in the outermost point of 
the convective core and the one in the innermost point of the 
radiative region. The final effect is to have a partial mixing 
at the edge of the convective core, and thus to remove the 
discontinuity in fj,. 

Furthermore, in models with a mass M ~ Ml cc , e.g. 
M — 1.2 Mq, the convective core is so small (m/M ~ 0.01) 
that the period spacing resembles the behaviour of the 1 
Mq model. We notice, however, the appearance of oscilla- 
tory components in AP in the model with X c ~ 0.1 (see 
Fig. I13[) . The sharp variation of N located at IIo/II r ~ 0.1 
is large enough to generate components with a periodicity 
of 10 k in AP. In more massive models, the fi gradient be- 
comes larger and so does the amplitude of the components 
in the period spacing (see e.g. Fig. 1 14f) . 

4-1.3 Models with a receding convective core. 

In models with shrinking convective cores the situation is 
much simpler. If M > M gc the dominant term in the be- 
haviour of the radiative gradient is the opacity and, since 
k oc (1 + X), V ra( j decreases with time as X decreases: the 
boundary of the convective core is displaced towards the 
centre. A receding convective core leaves behind a chemi- 
cal composition gradient that is responsible for an abrupt 
change in the N profile (as in Fig. [1} . Such a sharp feature 
which is a direct consequence of the evolution of a convec- 
tive core, leaves a clear signature in the periods of gravity 
modes. 

That behaviour is shown in Figures [K] and \W\ for mod- 
els of 1.6 and 6 Mo(for larger masses, the behavior is al- 
most identical to that of 6 Mq). The periodicity of the com- 
ponents in AP can be easily related to the profile of the 
Brunt- Vaisala frequency by means of expression (|14[) . For 
instance, the Afc = 7 in Fig. 1161 (lower panel, X c — 0.50), 
corresponds to the sharp signal in N at IIo/n r ~ 0.14. As 
the star evolves, the sharp feature in N is shifted to higher 
IIo/IIr, and when ITo/IT M ~ 0.5 a kind of beating occurs in 
the period spacing due to the fact that the sampling fre- 
quency is about half the frequency of the periodic compo- 
nent. 

The amplitude of the variation of AP as a function of 
the mode order is well reproduced by Eq. ((26} (compare 
e.g. Figs. 1161 and [8)l. but not by Eq. (|13[1 that predicts a 
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Figure 13. Behavior of the Brunt- Vaisala frequency (upper left 
panel), of the hydrogen abundance profile (upper right panel) 
and of the 1 = 1 g-mode period spacing in models of 1.2 Mq. 
We consider, as in all the following figures, several models on the 
main sequence with decreasing central hydrogen abundance (X c 
0.7, 0.5, 0.3 and 0.1). 



sinusoidal behaviour. However, having a simple analytical 
relation between the amplitude of the components and the 
sharpness of 8N/N is not straightforward from Eq. ((26}. 

It can also be noticed that oscillatory components of 
small amplitude occur already in zero age main-sequence 
stars with M > 6 Mq (see e.g. Fig. 1161 solid lines). Although 
a chemical composition gradient is not yet present in these 
models, the bump in the Brunt- Vaisala frequency due to an 
increase of the opacitjjf] at a temperature of ~ 3 x 10 6 K 
(IIo/n r ~ 0.8) is able to produce such a deviation from con- 
stant AP. It is not surprising that the effects of a sharp 
feature located near the surface can mi mic the effect of a 
pertur bation in the core: as shown by iMontgomerv et all 
( 2003]) the signature in high order g-modes of a perturbation 
in iV located at a normalized buoyancy radius tbv = IIo /n r 
is aliased to a signal whose source is located at 1 — rev- 
The signal shown in Fig \W\ could indicate a source at 0.2 
IIo/II r which is in fact approximatively an alias (1 — 0.2) 
of the source at 0.8 Ilo/n,.. The amplitude of this signal in- 
creases with the stellar mass as the contribution of this opac- 
ity bump becomes dominant in the behaviour of V ra d (and 
therefore of N). In fact for large enough stellar masses, a con- 



2 Mainly due to C, O, Ne and Fe tran sitions (see e.j 
iRogers fc Iglesiaslll992l ; ISeaton fe Badneljlioolh . 



Figure 14. Same as Fig |13l for 1.4 Mq models. 



vective shell can appear at a temperature ~ 3 x 10 K. The 
amplitude of such components is however less than 1000 s 
and therefore much smaller than the amplitude due to the 
chemical composition gradient at the edge of the convective 
core. 

The g-mode period spacing is clearly different depend- 
ing on the mass and age of the models. While in models with- 
out (or with a very small) convective core the mean value 
of the period spacing decreases with age (see Sec. l4.1.Tjl . for 
more massive models the period spacing does not signifi- 
cantly change with age. In stars with larger convective cores 
the chemical composition gradient is located at a larger frac- 
tional radius, giving a smaller contribution to J ^rdx', and 
thus to AP, as predicted from the asymptotic expression. 
In these models the age effect is made evident, however, 
through the appearance of the periodic signal in the period 
spacing, whose periodicity is directly linked to the chemical 
gradient left by the evolution of the convective core. 

4.2 Effects of extra-mixing 

The comparison between theoretical models and observa- 
tions clearly shows that the standard stellar modelling 
und erestimates the size of the centra l mix ed region (see 
e.g. lAndersen et al. | Il99d ; iRibas et al.l 120000 . This fact is 
generally accepted, but there is no consensus about the 
physical processes responsible for the required extra-mixing 
that is missing in the standard e volution models: over- 
I Schaller et al.lll992h . microscopic diffusion 



shooting (e.g. 



l|Michaud et all 



2004), rotationally induced mixing (see e.g. 
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Figure 15. Same as Fig |13l for 1.6 Mq models. 

iMaeder fc Mevnetllioool : iMathis et aDl2004l . and references 
therein), or mixing generated by propagation of internal 
waves (e.g. I Young fc Arnettll2005l ). The shape of the com- 
position transition zone is a matter of great importance as 
far as asteroseismology is concerned. In particular it signif- 
icantly affects the term appearing in the Brunt- Vaisala 
frequency and plays a critical role in the phenomenon of 
mode trapping. 

It is therefore evident that the size and evolution of the 
convective core, as well as the /j, gradients that it gener- 
ates, can be strongly affected by the occurrence of mixing 
processes. In the following paragraphs we study how these 
effects are reflected on the high order g-modes. We have 
computed models with overshooting, microscopic diffusion, 
turbulent mixing, and we have compared their adiabatic g- 
mode periods with those derived for models computed with- 
out mixing and with the same central hydrogen abundance. 

4-2.1 Overshooting 

Penetration of motions beyond the boundary of convective 
zones defined by the Schwarzschild stability criterion has 
been the subject of many studies in an astrophysical context 
(see e.g. lZahrj|l99lh . Unfortunately, features such as exten- 
sion, temperature gradient and efficiency of the mixing in the 
overshooting region cannot be derived from the local model 
of convection currently used in stellar evolution computa- 
tions. As a consequence, this region is usually described in a 
parametric way. In the models considered here the thickness 
of the overshooting layer ov is parameterized in terms of the 



Figure 16. Same as Fig |13l for 6.0 Mq models. In the homoge- 
neous ZAMS model (X <: ~ 0.7, full lines) a sharp variation in N 
located at IIo/n r ~ 0.8 generates a periodic signal of small am- 
plitude in the period spacing. Note that for this model, to make 
such a component more visible in AP, the y-axis scale has been 
modified. 



local pressure scale height H p : ov — (3 x (min(r cv , H p (r cv ))) 
(where r cv is the radius of the convective core and /3 is a 
free parameter). We assume instantaneous mixing both in 
convective and in overshooting regions. The temperature 
gradient in the overshooting region is left unchanged (i.e. 
V = V ra d). Therefore overshooting simply extends the re- 
gion assumed to be fully mixed by convection. The larger 
hydrogen reservoir, due to an increase of the mixed region, 
translates into a longer core-hydrogen burning phase. 

The adopted amount of overshooting also determines 
the lowest stellar mass where a convective core appears. For 
sufficiently large values of /3, the convective core that devel- 
ops in the pre-main sequence phase persists during the main 
sequence in models with M <Ml cc (as in Fig. |9}. In these 
models the convective core is maintained thanks to the con- 
tinuous supply of 3 He that sustains the highly temperature- 
dependent nuclear reaction 3 He + 3 He — > 4 He + 1 H + 1 H, 
keeping the pp chain in an out-of-equilibrium regime. The 
inclusion of overshooting changes the value of the mass cor- 
responding to the transition between models with a convec- 
tive core that grows/shrinks during the main sequence (Fig. 
[9] right panel). And finally, the effect of overshooting on the 
H gradients depends on whether the nuclear reactions occur 
only inside the convective core or also outside. 

In Figure [17] we present the chemical composition pro- 
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Figure 17. Behavior of the Brunt- Vaisala frequency (upper row), of the hydrogen abundance profile (central row) and of the 1 = 1 
g-mode period spacing (lower row) in models of 1.2,1.6 and 6 Mq computed with (thick lines) or without (thin lines) overshooting. 



file, the behaviour of TV and of the period spacing, in models 
computed with overshooting (f3 = 0.2). These models have 
a larger fully mixed region than those computed without 
overshooting. The chemical composition gradient is then dis- 
placed to a higher mass fraction. If we compare with models 
of similar central hydrogen abundance, however, this does 
not necessarily imply that the sharp feature in N is located 
at a different normalized buoyancy radius (IIo/II^). 

• In 6 M© models (right column in Fig. I17[l , for instance, 
neither the sharpness of the abrupt variation in TV, nor its 



location in terms of II M , change when comparing models 
computed with and without overshooting. 

• The situation changes in lower mass models, e.g. in 
1.6 Mq models (central column in Fig. I17|l . Here the pe- 
riodicity of the components in AP differs if we include over- 
shooting or not. A change in the location, but also in the 
value of the p gradient (as nuclear reactions take place out- 
side the core as well) , is responsible for a different behaviour 
of the oscillatory components in AP. 

• In models with a mass M ~ Ml cc , e.g. M = 1.2 Mq, 
the inclusion of overshooting dramatically increases the size 
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Figure 18. Behavior of the Brunt- Vaisala frequency in models of 
1.6 Mq with X c ~ 0.3. The different lines correspond to models 
calculated with no extra-mixing (continuous lines), overshooting 
(dashed) and helium diffusion (dotted). The different location 
and sharpness of the chemical composition profile determine the 
behaviour of N. 



of the fully mixed region. The oscillatory components in AP 
have different periods and much larger amplitudes than in 
the case of models without overshooting (see left column of 

Fig. ni]i. 

As we mentioned above, not only the extension of the 
overshooting region is uncertain, but also its temperature 
stratification. However, the effect on AP of considering con- 
vective penetration (V = V a d in the "extended" convective 
core) instead of simple overshoot ing is found to be small (see 
IStraka et al. l l2005l ; lGodartll2007l ). 



4-2.2 Effects of microscopic diffusion 

Other physical processes, different from overshooting, can 
lead to an increase of the central mixed region or modify 
the chemical composition profil e near the core . 

iMichaud et""ai1 l|2004l ) and iRichardl (|2005l ) have shown 
that microscopic diffusion can induce an increase of the con- 
vective core mass for a narrow range of masses, from 1.1 to 
1.5 Mq, and that the effect decreases rapidly with increas- 
ing stellar mass. In this mass range, as previously described, 
the mass of the convective core increases during the MS evo- 
lution instead of decreasing, as it occurs for larger masses. 
As a consequence, a sharp gradient of chemical composition 
appears at the border of the convective core, making the 
diffusion process much more efficient in that region. 

Figure[TH]shows the effect of including gravitational set- 
tling, thermal and chemical diffusion of hydrogen and helium 
in a 1.6 Mq model. If only He and H diffusion is considered in 
the modelling we find that a smoother chemical composition 
profile is built up at the edge of the convective core (see Fig. 
I18[) preventing the presence of a discontinuity in /i caused 
by a growing convective core. Comparing models calculated 
with or without diffusion, the location of the convective core 
boundary is not significantly changed. Nevertheless, a less 



Figure 19. Period spacing of I = 1 g modes as a function of the 
radial order k for the 1.6 Mq models presented in Fig. 1181 



sharp variation in the Brunt- Vaisala frequency is respon- 
sible for a considerable reduction of the amplitude of the 
components in the period spacing (see Fig. 1190 . It was in 
fact predicted in Section [3.11 that a discontinuity not in TV 
itself, but in its first derivative, would generate a periodic 
component whose amplitude decreases with the period: this 
simplified approach is therefore sufficient to account for the 
behaviour of AP in models with a smoother chemical com- 
position gradient (see Eq. I16|l . 

When diffusion of Z is also included, the overall effect 
of diffusion is no longer able to erode the sharp chemical 
composition gradient and to prevent the formation of a semi- 
convective zone, on the contrary, diffusion of Z outside the 
core makes the occu r rence of semi-convect i on easier (see also 
Richard etail l200ll ; IMichaud et~ail |2004| ; iMontalban et al.l 



20071) . 



If chemical diffusion is accountable for a smoother 
chemical composition profile in intermediate- and low-mass 
stars, we find that such an effect disappears as higher masses 
are considered and evolutionary time-scales decrease. In fact 
we find that in models with M > 4 Mq diffusion has no 
effect on the Brunt- Vaisala frequency profile near the hy- 
drogen burning core, nor on the behaviour of the period 
spacing. 

As higher masses are considered, the effect of micro- 
scopic diffusion on the stellar structure near the core be- 
comes negligible but, other mixing processes can partly 
erode the chemical composition gradients. 



4-2.3 Rotationally induced mixing 

Rotationally induced mixing can influence the internal dis- 
tribution of n near the energy generating core. Different 
approaches have been proposed to treat the effects of ro- 
tation on th e transport of angular momentum and chemi- 
cals (see e.g.lMaeder fc Mevnetboool ; iHeger fc Langerll2000l ; 
iPinsonneault et al.lll989l and references therein). Such an 
additional mixing has an effect on the evolutionary tracks 
which is quite similar to that of overshooting, but it leads 
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Figure 20. HR diagram showing evolutionary tracks of 6 
Mq models calculated with different extra-mixing processes. The 
evolutionary track computed with Dt = 500 is superposed to 
that without extra-mixing (continuous line). 



also to a smoother chemical composition profile at the edge 
of the convective core. 

Since our stellar evolution code does not include a con- 
sistent treatment of rotational effects, we simply include the 
chemical turbulent mixing by adding a turbulent diffusion 
coefficient (Dt) in the diffusion equation. In our parametric 
approach Dt is assumed to be constant inside the star and 
independent of age. 

The simplified parametric treatment of rotationally in- 
duced mixing used in this work has the aim of showing that 
if an extra-mixing process, different from overshooting, is 
acting near the core it will produce a different chemical com- 
position profile in the central regions of the star and leave 
a different signature in the periods of gravity modes. The 
effects of such a mixing on the evolutionary tracks (see Fig. 
I20[) and on the internal structure clearly depend on the value 
of D T . 

As an example we consider models of a 6 Mq star com- 
puted with three different values of the turbulent-diffusion 
coefficient D T = 5 x 10 2 , 5 x 10 3 and 5 x 10 4 cm 2 s" 1 . The 
lowest Dt value provides an evolutionary track that over- 
laps the one computed with no mixing, however, its effect on 
the chemical composition gradient (see Fig. I2ip is sufficient 
to affect AP for high order modes. In order to significantly 
change the periods of low-order modes a much more effective 
mixing is needed. The value Dt = 5 x 10 3 leads to a slightly 
more luminous evolutionary track, but such a mixing has 
a substantial effect on the period spacing: the amplitude of 
the periodic components in AP becomes a decreasing func- 
tion of the radial order k (see Figures [21] and I22|l . As in 
the case of helium diffusion (see Sec. I4.2.2[) this behaviour 
can be easily explained by the analytical approximation pre- 
sented in Sec. 13.11 (Eg. I16|l . provided the sharp feature in N 
is modelled not as a step function but, for instance, as a 
ramp (Eg. [To)) . 

If a significantly more effective mixing is considered (e.g. 
Dt = 50000, see Fig. 1201 and I23[) , the corresponding evolu- 
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Figure 21. Behavior of the Brunt- Vaisala frequency (upper left 
panel), of the hydrogen abundance profile (upper right panel) and 
of the 1=1 g-mode period spacing in models of 6.0 Mq computed 
with a turbulent diffusion coefficient Dt = 500. 



tionary track is close to that obtained by including a classi- 
cal overshooting but the periodic components in AP are no 
longer present. 

The effects of such a turbulent mixing on low-order 
gravity modes and avoided crossings will be addressed in 
detail in a future work. 

In order to check that our parametric approach is at 
least in qualitative agreement with the outcome of models 
where rotationally induced mixing is treated consistently, 
we present for comparison a seque nce of models computed 
with the Geneva evolutio nary code (jMevnet fc Maedenl2000l : 
lEggenberger et al]|200^ . We considered a 6 Mq model with 
an initial surface rotatio nal velocity of 25 km s _1 , the typi- 
cal vsini for SPBs stars (jBriauet et alJl2007t ) . As the model 
evolves on the MS, the effects on the HR diagram (Fig. I24|l 
and on the chemical composition gradient in the core (Fig 
I25p are very similar to those of the model computed with a 
uniform turbulent diffusion coefficient Dt = 5 x 10 3 . This 
is not surprising since, in the central regions, the total tur- 
bulent diffusion coefficient shown in Fig. [55] does not change 
considerably in time and its magnitude is of the order of a 
few thousands cm 2 s _1 . 



4.3 Modes of different degree and low radial order 

In the previous sections we analyzed the properties of AP 
for 1 = 1, high order g-modes. Does the period spacing com- 
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Figure 22. As in Fig. [22] but for models with D T = 5000 



Figure 23. As in Fig. |22]but for models with D T = 50000. The 
effect of such a mixing on the evolutionary track is shown in Fig. 



puted with modes of different degree £, or with low radial 
order k, have a similar behaviour? In Fig. [27] we consider 
the period spacing computed with t = 2 periods Pk,i=2 
scaled by the ^-dependent factor suggested by Eq. [4] i.e. 
Pfc 1-2 = v3 Pk,t=2- Except for the oscillation modes with 
the lowest radial orders, we notice in Fig.[27]that the period 
spacing of modes of degree I = 1 and 2 have the same be- 
haviour, provided that the dependence on I given by Eq. [4] 
is removed. 

Though the asymptotic approximation is valid for high 
order modes, we see in Fig. [27] (as well as in the figures pre- 
sented in the previous sections) that the description of the 
period spacing as the superposition of the constant term de- 
rived from Eq. [4] and periodic components related to V M , is 
able to accurately describe the periods of gravity modes even 
for low orders k. This suggest that, at least qualitatively, the 
description of g-modes presented in this work could repre- 
sent a useful tool to interpret the behaviour of low order 
g-modes observed in other classes of pulsators, such as ft 
Cephei and S Scuti stars. This subject will be addressed in 
a second paper. 



5 CHALLENGES FOR ASTEROSEISMOLOGY 

Asteroseismology of high-order g-mode main-sequence pul- 
sators is not an easy task. The long oscillation periods and 
the dense frequency spectrum in these stars require long and 
continuous observations in order to resol ve single oscillation 
frequencies. As clearly stated already in iDziembowski et al.l 
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Figure 24. HR diagram showing evolutionary tracks of 6 
Mq models calculated with the Geneva evolutionary code. The 
full line evolutionary track is computed without rotation, whereas 
in the models evolving on the dotted track an initial surface ve- 
locity of 25 km s _1 is assumed. 



Probing the properties of convective cores through g modes 15 



O.E 



0.6 



x 0.4 



0.2 



I 1 1 1 1 I 



6 M„ 



/ 



V, = 25 km/s- 

0.1 0.2 0.3 0.4 0.5 
m/M 
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with an initial surface velocity of 25 km s . 
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Figure 26. Total diffusion coefficient Dx as a function of the nor- 
malized mass in 6 Mq models evolving on the main sequence. The 
initial surface velocity assumed is 25 km s _1 . Dx includes both 
the effects of shear (D^a,.) and meridional circulation (D c g). 



l|l993h , observations of at least 70 days are needed to resolve 
the period spacing in a typical SPB star. In addition to large 
observati onal efforts that have been made from th e ground 
(see e.g. |Pe Cat fc Aertsl l2002l . |De Cat et afll2007l for SPB 
stars and Mathias et alj 12004 and Henry et a!T 2007 for 7 
Dor), long, uninterrupted space-based photometric time se- 
ries will soon be provided by C0R0T. In order to estimate if 
the C0R0T 150 day-long observations have a frequency res- 
olution sufficient to resolve the period spacing of a typical 
SPB star, we simulated a time series of I — 1 and £ — 2 
high-order g modes that are expected to be excited in a 



Figure 27. Period spacing for a 6 Mq models with Xc = 0.5. 
The period spacing computed with 1 = 2 modes is multiplied by 
\/3, as suggested by Eq. [4] 



6 Mq star with Xc = 0.2 computed without overshooting. 
The excitation of oscillation m odes has been comp uted using 
the non-adiabatic code MAD IIDupret et alj .2003). Random 
initial phases have been considered in the sine curves de- 
scribing the 25 frequencies included in the simulations, and 
amplitudes of 1 and 0.5 (on an arbitrary scale) have been 
assumed for £ = 1 and 1 = 2 modes respectively. The gen- 
erated time series has been analyzed using the package Pe- 
riod04 (|Lenz fc Bregerll2005h . As it is shown in the resulting 
power spectrum (see Fig. [28}, for such long observational 
runs (150 days) the input frequencies can be accurately re- 
covered even for the largest periods of oscillation. Thanks 
to the high frequency resolution, the departures from a con- 
stant period spacing are also evident in the power spectrum. 

We then simulate a time series for an additional 6 
Mq model that, despite having the same surface properties 
(log(L/Lo) = 3.28, log (T c ft) = 4.22) as the previous one, 
is computed with turbulent diffusion (Dx = 5000). In this 
case the power spectrum (see Fig. I29|l shows a regular period 
spacing that, on the basis of the high frequency resolution, 
can be easily distinguished from the one described in Fig. 
[281 

Even though frequency resolution may no longer be 
an issue in the very near future, there is a second well 
known factor limiting asteroseismology of SPB and 7 Dor 
stars: the effects of rotation on the oscillation frequency can 
severely complicate the high-or der g-mode spectra . Refer- 
ring once more to the work by iDziembowski et al.1 l|l993h . 
a first requirement in order to treat rotational effects as 
perturbations on the oscillation frequencies is the angular 
rotational velocity being sufficiently smaller than the oscil- 
lation frequency (v ro t <§C 2nR/P, where R is the radius of 
the star and P the oscillation period). In the case of the 
models of an SPB star considered previously in this sec- 
tion, this translates into v ro t <siC 100 km s -1 considering 
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Figure 28. Power spectrum of simulated time series for a 6 
Mq model. The vertical solid and dashed lines represent the input 
oscillation frequencies of, respectively, 1=1 and 1 = 2 modes. 



the modes of longest period: this is significantly larger than 
the avera ge vsini (~ 25 km s~ ) measured in SPB stars 
(see e.g. iBriquet et al][2 007). In the case of 7 Dor stars 
(see e.g. ISuarez et alj|200a ) a similar estimate limits the va- 
lidity of the perturbative approach to rotational velocities 
of ~ 50 — 70 km s _1 : for fast er rotators non-perturbative 
approaches are needed (s ee iDintrans fc Rieutordl l2000l: 
iRieutord fc Dintransj[2002l ). 

Even in the case of slow rotators, however, the rota- 
tional splittings may become as large as the period spac- 
ing itself. Such a large effect of rotation on the frequency 
spectrum can severely complicate the identification of the 
azimuthal order m and the degree £ of the observed modes. 
In fact, a simple mode identification based on the regular 
pattern expected from high-order g-modes becomes inappli- 
cable if the rotational s plitting is of the sa me order as AP 
(see e.g. Fig. 13 in lDziembowski et al.lll993T ). The identifica- 
tion of the modes would then need to be provided by pho- 
tometric and spectroscopic mode identification techniques 
(see e.g. lBalonalll986l:lGarrido et al.lll990l;lAerts et alJ |l992: 



Mantegazzdl200d ; lBriauet fc Aertj2003l ; lDupret et alJhoOS 
Zimall2006l ). 



6 SUMMARY AND CONCLUSIONS 

In this work we investigated in detail the properties of high- 
order gravity in models of main-sequence stars. The chemical 
composition gradient that develops near the outer edge of 
the convective core leads to a sharp variation of the Brunt - 
Vaisala, frequency. As we presented in Sec. [2] the latter is 
responsible for a periodic trapping of gravity modes in the 
region of chemical composition gradient, and it directly af- 
fects the period spacing of g modes. 

In analog y with the works on white dwarf s by 
iBrassard et al.1 (Il992l ) and iMontgomerv et aL I i|2003h . we 
show that in the case of main sequence stars analytical ap- 
proximations can be used to directly relate the deviations 
from a uniform period spacing to the detailed properties of 
the /i-gradient region that develops near the energy gener- 
ating core. We find that a simple approximation of g-mode 




Figure 29. Same as Fig. 1281 for a 6 Mq model computed with 
turbulent diffusion in the core. 



periods, based on the variational principle of stellar oscil- 
lations, is sufficient to explain the appearance of sinusoidal 
components in the period spacing. This approximation (see 
Sec. 13. ip relates the periodicity of the components to the 
normalized buoyancy radius of the glitch in N, and the am- 
plitude of the components to the sharpness of the feature 
in N. In particular, if the sharp variation in N is mod- 
elled as a step function, the amplitude of such components 
is expected to be independent of the order of the mode k; 
whereas if the glitch in N is described with a ramp, the 
amplitude of the components decreases with k. A more ac- 
curate semi-analytical approximation of the period spacing, 
which considers the effects of the sharp feature in N on the 
eigenfunctions, is also given in Sec. 13.21 

We then presented a survey of the properties of high- 
order g modes in main sequence models of masses between 1 
and 10 Mq (see Sec.|4]). As a general result we found that, in 
models with a convective core, the period spacing of high- 
order g modes is accurately described by oscillatory com- 
ponents of constant amplitude, superposed to the mean pe- 
riod s pacing predicted by the asymptotic theory of iTassoull 
(1980). In Sec. 14.11 we showed that the period spacing de- 
pends primarily on the extension and behaviour of the con- 
vective core during the main sequence and, therefore, on the 
mass of the star. 

In models without a convective core (see Sec. I4.1.ip the 
mean AP considerably decreases during the MS, whereas 
no significant deviation from a constant period spacing is 
present. For an intermediate range of masses (see Sec. 14. 1.2)l 
the convective core grows during most of the MS, gener- 
ating an "unphysical" discontinuity in fj, if no mixing is 
added in the small semiconvective region that develops. 
We find that the behaviour of AP, and in particular the 
appearance of periodic components, depends on the treat- 
ment of this region. It is interesting to notice that 7 Do- 
radus stars are in the mass domain where models show a 
transition between growing to shrinking convective cores on 
the main sequence. Gravity modes could therefore repre- 
sent a valuable observational test to discriminate between 
the different prescriptions used in stellar models (see e.g. 
iPopielski fc Dziembo wski 2005) to introduce the required 
mixing at the boundary of the convective core. In models 
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with higher masses, the convective core recedes during the 
main sequence (see Sec. I4.1.3[) : this leaves behind a jj, gra- 
dient that generates clear periodic components in AP. We 
found that the analytical expression derived in Sec. 13,11 al- 
lows to accurately recover the location and sharpness of the 
H gradient from the amplitude and periodicity of the com- 
ponents in AP. In this mass domain, though the average 
period spacing does not change substantially with the age, 
the periodicity of the components does, and it therefore rep- 
resents an indicator of the evolutionary state of the star. 

In Sec. 14.21 we showed that also extra-mixing processes 
can alter the behaviour of AP, since they affect the size 
and evolution of the convective core, as well as the sharp- 
ness of the \i gradient. We first compared models with the 
same X c , but computed with and without overshooting (see 
Sec. I4.2.T1I . We found that in models with small convective 
cores, or where nuclear reactions take place also in the ra- 
diative region, the different size of the fully-mixed region 
changes the periodicity of the components in AP. In Sec. 
14.2.31 we described how chemical mixing can severely affect 
the amplitude of the periodic components in AP. In models 
where turbulent mixing induced by rotation is considered, 
the smoother /i profile near the core leads to a discontinuity 
not in A^ itself, but in its first derivative: as suggested by the 
analytical approximation in Sec. 13.11 this leads to periodic 
components in AP whose amplitude decreases with the or- 
der of the mode. In the case of SPB stars, in particular, we 
find that the mixing induced by the typical rotation rates 
observed (i.e. ~ 25 km s _1 ), is sufficient to alter significantly 
the properties of the g-mode spectrum. 

Finally in Sec [S] we discussed the difficulties encountered 
in the asteroseismology of 7 Doradus and SPB stars. Even 
though a frequency resolution sufficient to resolve closely 
spaced periods will be provided by the forthcoming space- 
based observations, an asteroseismic inference on the inter- 
nal structure will only be possible for stars with very slow 
rotation rates, and with reliably identified pulsation modes. 
Once these conditions are reached, we will be able to ac- 
cess the wealth of information on internal mixing which, as 
shown in this work, is carried by the periods of high-order 
gravity modes in main-sequence objects. 
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